Origin of combination frequencies in quantum magnetization oscillations of 

two-dimensional multiband metals 



Thierry Champel 
Commissariat a VEnergie Atomique, DRFMC/SPSMS 
17 rue des Martyrs, 38054 Grenoble Cedex 9, France 

(February 1, 2008) 



The influence of chemical potential oscillations on the magnetization oscillations in two- 
dimensional multiband metals is investigated. In a first part, the analytical derivation of Alexan- 
drov and Bratkovsky [Phys. Rev. B 63, 033105 (2001)] of magnetic quantum oscillations in two- 
dimensional multiband metals with a fixed number of electrons is shown to be mathematically 
incorrect ; the chemical potential oscillations appearing in the arguments of the Fourier components 
were not taken into account. In a second part, we derive an approximative Fourier series of the 
magnetization oscillations in the regime of small chemical potential oscillations. The main result is 
that combination frequencies with significant amplitudes are found if the individual band frequencies 
differ significantly. 



The observation of quantum magnetization oscillations (the de Haas- van Alphen effect) is one of many manifestations 
of the quantization of the electronic energy spectrum in a magnetic field H. This quantum effect is observable both in 
two-dimensional (2D) and three-dimensional (3D) metals. However, on account of the importance of chemical potential 
oscillations, the de Haas- van Alphen (dHvA) effect in 2D metals is quite different from its 3D counterpart. It becomes 
dependent on experimental conditions : the magnetization oscillations show different characteristics corresponding to 
either imposing a constant number of electrons (for which the chemical potential /i oscillates with the magnetic field) 
or to imposing a fixed chemical potential 1-3 . For example, in 2D multiband metals, unusual combination frequencies 
of the ordinary frequencies of the individual bands appear in the Fourier spectrum of magnetization oscillations in 
the presence of chemical potential oscillations. The observation 4 of these additional contributions to the dHvA effect 
with significant amplitudes enables a quantitative comparison to theory to be made. For this purpose, an analytical 
description of the dHvA effect in 2D multiband metals including chemical potential oscillations on a quantitative basis 
is needed. 

Recently, Alexandrov and Bratkovsky 5 have proposed an analytical derivation for the amplitudes of the combination 
frequencies in magnetic quantum oscillations. They argued that, for a fixed number of electrons TV, the combination 
frequencies originate from the oscillation of a squared term, only present in the expression of the free energy F in 
the canonical ensemble. After calculating the grand canonical thermodynamic potential f2, they obtained the explicit 
form of the free energy [Eq. (13) of Ref. 5] using the relation F = fl + fj,N. From their expression (13), they got 
straightforwardly the Fourier harmonics of the combination frequencies. 

According to Ref. 5, the field oscillating part F of the free energy F is related to the field oscillating part Cl of the 
grand canonical thermodynamic potential f2 through 



where p is the total density of states. Following Alexandrov and Bratkovsky 5 , it is the second nonlinear oscillatory term 
in Eq. (1), specific to the canonical ensemble, which produces the combination frequencies. However, it is important 
to stress here that il is an explicit function of the variable fi with the ratio /J,/to ca appearing in the arguments of the 
Fourier components (where uj ca is the cyclotron frequency of the individual band a), see explicit form Eq. (6), Ref. 
5. For a fixed total number of electrons N, the correct treatment is thus to eliminate the variable fi in the expression 
for SI via the condition 



Under a magnetic field, the chemical potential consists of a constant part fio (the value at zero magnetic field) plus 
an oscillating part /2 given by the implicit equation 
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In 3D metals, the chemical potential oscillations are negligibly small 6 since p ~ uj c \] u) c / po and po 3> ^ c (in this 
part, the band index a is omitted for convenience). The oscillating part p appears in the arguments of the Fourier 
components of fl in the combination p/cj c ~ ^/cj c /po- Therefore, the direct substitution of the variable p by the 
constant part po in the explicit expression of Ct is valid 6 : Sl(/Lt) w f2(/xo). It means that the exact resolution of implicit 
Eq. (3) is not needed. In this case, the explicit expression for can effectively be considered as a Fourier series . 
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Moreover, the second term of Eq. (1), ydSl/dpj /2p ~ pLo 2 c (u c / po) , is negligible compared with il ~ puj 2 ^uj c /pQ. 

In 3D multiband metals, Eq. (1) is therefore applicable but combinations frequencies arising from the squared term 
in expression (1) have very small amplitudes. 

In 2D multiband metals and at very low temperature T <C u) c , the situation is different. As expected in these 
conditions, the amplitude of chemical potential oscillations is much greater than in 3D metals. According to Eq. 

(13) of Ref. 5, ~ pcu 2 ~ (j)Q,/dpj /2p at low temperature and for a weak impurity scattering : the two terms 

in expression (1) are then of the same order. However, at the same time, we have p ~ lo c so that p/uj c = 0(1). 
Consequently, the direct substitution of the variable p/oj c by Pq/lu c in the sine arguments of the explicit expression 
for is no longer valid. Contrary to the 3D metals, the exact resolution of the implicit Eq. (2) is needed. The 
consideration of the explicit expression for F as a Fourier series, as done in Ref. 5, is now incorrect. Furthermore, on 
account of the non negligible oscillating part p, the individual band contributions are also mixed via the nonlinear 
coupling in f2, so that combination frequencies can be produced by both terms in Eq. (1). Hence, we conclude that 
the analysis proposed by Alexandrov and Bratkovsky 5 for explaining the appearance of combination frequencies in 
2D multiband metals is mathematically incorrect. 

Moreover, the mechanism proposed by Alexandrov and Bratkovsky is based on the relation F = SI + pN which 
is valid only at the thermodynamic limit. In this limit, the thermodynamic quantities evaluated in the canonical 
ensemble or in the grand canonical ensemble are equal for the same given experimental condition. The magnetization 
oscillations function M is then equivalently derived by 

M = - 

Consequently, the mechanism responsible for the combination frequencies in the Fourier spectrum of magnetization 
oscillations can not a priori depend on the way the magnetization is derived, that is to say on the use of a specific 
thermodynamic potential. 

Our following goal is then to point out how the combination frequencies arise by considering directly the expression 
for the magnetization oscillations in the relevant thermodynamic limit. Exactly like in the 3D Lifshitz-Kosevich deriva- 
tion 6 , the calculations of magnetization oscillations 8 are technically easier with the grand canonical thermodynamic 
potential Q. Therefore, the magnetization oscillations are derived by way of 
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which is a function of the grand canonical variables p, V and T. For a fixed total number of electrons N, the difficulty 
is still to eliminate the variable p in the expression for M via the condition (3). As previously stated, the expression 
(3) for the chemical potential oscillations is the key-equation to understand the difference between 2D and 3D dHvA 
effects. Contrary to the 3D case, its exact resolution is needed. The possible presence of a finite intrinsic reservoir 
of electrons 2 which reduces the chemical potential oscillations has also to be taken into account as an additional 
parameter. As a result, in 2D metals, the magnetization oscillations become more or less sensitive to the presence of 
the chemical potential oscillations and exhibit significantly different behaviors depending on the presence or not of 
a finite reservoir of electrons 2 ' 9 ' 10 . In 2D multiband metals this higher sensitivity to chemical potential oscillations 
is expressed by the presence of the combination frequencies in the Fourier spectrum of magnetization oscillations at 
very low temperatures, as shown numerically by Nakano 11 . Here, our aim is to prove analytically the existence of 
these combination frequencies in the magnetization oscillations. 

In 2D one-band metals, the oscillating part p and the magnetization oscillations M are linked through a simple 
relation of proportionality so that equation (3) and equation (4) can be compiled to give a single equation 9 
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where lu c = eH/mc is the cyclotron frequency and m the effective mass, Mq is the magnetization at saturation, 
A, = 2tt 5 'Ik bT/HlOc, T is the width of the Landau levels due to impurity scattering (assuming a Lorentzian broadening) , 
and p e is the electron's magnetic moment. The dimensionless parameter R measures the strengh of the coupling to 
the background reservoir 2 (for a more precise definition of R, see Ref. 9). The chemical potential oscillations are 
responsible for the nonlinearity of Eq. (5). For R » 1, they are strongly reduced : then, formula (5) yields directly 
the Fourier development of the magnetization oscillations in 2D metals for a fixed chemical potential 1 . 
In the multiband case, Eq. (3) and (4) give for ksT <C p and hoj c <C p 
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where p a = p Qa + p, p 0a = p — A a0 is independent of the magnetic field, M Qa = p a po, 
field density of states per spin projection in the band a with the effective mass m a , cj ca — eH/m. 
R = pr/ J2 a Pa (with PR the density of states of a non-quantized background reservoir), and 
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We note here that the equations (6) and (7) can not be combined to give a single equation as was the case for 
one band. Nevertheless, the effects of p on the magnetization oscillations are qualitatively the same. As for the 
one-band metal, for R ^> 1, the chemical potential oscillations are strongly damped and the Fourier expression with 
classical frequencies f a = 2ircm a po a /e is recovered. For a smaller value of R, the contributions of each band to the 
magnetization oscillations are mixed through the oscillating part of the chemical potential p. However, from Eq. (6) 
and (7), it is not clear whether a Fourier analysis for the magnetization oscillations is relevant. Let's try to make a 
Fourier development. For this purpose, in Eq. (6) and (7) we separate the two different parts in the sine arguments : 
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At any finite temperature or impurity scattering, the quantity p/Tiio ca is strictly less than 0.5(1 + R)~ 1 m a /m and is 
reduced further at higher temperature or lower magnetic field. Thus, expanding p/fkj ca in powers of a temperature or 
impurity reduction factor, Eq. (9) could be solved by iteration. However, the nonlinearity makes the resolution of this 
self-consistent equation somewhat cumbersome, especially for strong chemical potential oscillations \p/hco ca \ ~ 1/2. 
According to Eq. (8) and (9), it is worth noting that, generally, for a fixed number of electrons, a Fourier development 
of the magnetization oscillations may not exist. Nevertheless, in some particular regimes or under specific conditions, 
the description in terms of a Fourier series seems possible locally, that is to say for a finite range of magnetic field 
(which depends on temperature, on the parameter i?, on the impurity broadening T and also on the electron hopping 
integral t in quasi-2D metals). Indeed, in the regime of small but non negligible chemical potential oscillations 
\p/huj ca \ < 1/27T, the linearization of Eq. (8) and (9) is possible (for the first significant /) and yields at first order in 
chemical potential oscillations 
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The substitution of (11) into (10) leads to the Fourier series expansion 
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In the presence of several bands, the assumption of small chemical potential oscillations is not so restrictive even for 
R = 0, since the amplitudes of oscillations are reduced by the extra factor m Q /m < 1. In this regime, it is the second 
term of (12) which is responsible for the presence of the combination frequencies / = lf a ± l'f a >. Their amplitudes 
are 

^kfl^M 0a ±l'^M 0a ,). (13) 
2(1 + R) \ m to J 

At T — 0, the ratio of the amplitude of the combination frequency f a + f a ' harmonic to the single band a amplitude 
given by the first term is (A\ = 2/ir without spin-splitting and impurity factors) 

1 m a f 1+ U\ 



1 + R m \ f a 

This value is not incompatible with the condition of linearization. It shows that the presence of combination frequencies 
is significant if f a and f a > do not have the same order of magnitude. In experiments of Shepherd et al. 4 , the ratio of 
band frequencies is of the order of 4, which is consistent with this condition for the observation of the combination 
frequencies. In the numerical work of Nakano 11 , significant amplitudes for the combination frequencies are found 
in the same configuration of f a i / f a ~ 4. It means that the 2D multiband metals with individual band frequencies 
such that f a > lfa 3> 1 are ideal systems to observe combination frequencies in the Fourier spectrum of magnetization 
oscillations at very low temperatures. 

Moreover, from Eq. (12) it is seen that the chemical potential oscillations modify the amplitudes of the single 
band frequencies. Indeed, the amplitude of the frequency / — Lf a consists of the usual amplitude (first term in the 
right-hand side of Eq. (12)) plus an infinity of terms produced by the chemical potential oscillations in the second 
term when L = I ± /' . For example, the ratio of the amplitude of the harmonic L (frequency Lf a ) in the presence of 
small chemical potential oscillations to the usual one is 
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For the first harmonic L = 1, this ratio becomes 
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Using the formula 1/1(1 + 1) = 3/2, we find at zero temperature and in absence of spin-splitting and impurity 
scattering : 

3^TO^ 
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In the regime of small chemical potential oscillations, the ratio (17) deviates slightly from unity, so that the ratio of 
the amplitude of the combination frequency f a + f a < to the single band f a amplitude [Eq. (14)] is unaffected by this 
correction factor. However, the expressions (15), (16) and (17) derived with the linearisation of the chemical potential 
oscillations indicate that in the regime of strong chemical potential oscillations (for which numerical calculations of 
the magnetization oscillations are unavoidable) the usual single band amplitudes are also significantly affected 2 ' 10 . 

In conclusion, the dHvA effect in 2D multiband metals at a constant number of electrons has been investigated 
analytically, emphasizing the clear difference from 3D metals. In the presence of strong chemical potential oscillations, 
the quantitative description of the magnetization oscillations in terms of a Fourier series may break down and has 
to be done numerically. In the regime of small chemical potential oscillations, the Fourier analysis is still a good 
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approximation. In 2D multiband metals, these oscillations are responsible for the presence of combination frequencies 
which can occur with a significant amplitude if the individual band frequencies differ significantly. 
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